NACA TN 2245 


CASE FILE 
COPY 

NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 2245 

CALCULATION OF COMPRESSIBLE POTENTIAL FLOW PAST SLENDER 

BODIES OF REVOLUTION BY AN INTEGRAL METHOD 

By Milton M. Klein and W. Perl 

Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 



Washington 

December 1950 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2245 


CALCULATION OF COMPRESSIBLE POTENTIAL FLOW PAST S LEN DER 
BODIES OF REVOLUTION. BY AN INTEGRAL METHOD 
By Milton M. Klein and W. Perl 


SUMMARY 

An integral method, previously used to obtain compressible flow 
past two-dimensional shapes, has been applied to the calculation of 
compressible potential flow past slender bodies of revolution. Good 
agreement of the resulting velocities on ellipsoids of revolution 
with- those obtained by other methods was found. Flow -field conditions', 
such as the extent of local supersonic regions, were calculated and 
discussed. The free-stream Mach number beyond which a continuous 
potential -flow solution ceases to exist was found to be very close to 1. 


INTRODUCTION 

The partial differential equation for the compressible potential 
flow past a body of revolution is considerably more difficult to solve 
than the corresponding plane (two-dimensional) equation. As a conse- 
quence, fewer extensions of low-speed or incompressible -flow solutions 
to higher subsonic speeds have been made for bodies of revolution than 
for two-dimensional shapes. The principal work on bodies of revolution 
seems to have been confined to solutions of the partial differential 
equation after it has been linearized under the assumptions of small 
perturbations and of free-stream Mach number not too close to 1. For 
example, a treatment of subsonic flow by the method of sources and 
sinks is given in reference 1. The case of wholly supersonic flow,, 
which has been more extensively treated, is discussed in reference 2. 
These solutions may be conveniently expressed in the form of* a 
similarity rule ( Brandt 1-Glauert rule for bodies of revolution) that 
gives the flow pattern for a given Mach number if the flow pattern 
at another Mach number is known. Correct forms of the rule are pre- 
sented in references 3 to 5. 

The solution of the linearized differential equation yielding 
the Prandtl-Glauert rule may be regarded, as in the two-dimensional 
case, as the first term in an expansion of the potential in powers of 
a thickness -ratio parameter. The amount of labor required to obtain 
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higher -order terms of the expansion, is considerable for a two- 
dimensional body and much greater for the corresponding body of 
revolution. Only one calculation, for the ellipsoid of revolution, 
of the second term of such an expansion appears to have been made 
(reference 6) . ' 

Because of the large amount of work and the mathematical diffi- 
culties involved in obtaining higher approximations for bodies of 
revolution through the differential -equation approach, it is particu- 
larly desirable to survey other possible approaches to the problem of 
calculating the flow past a body of revolution at high subsonic Mach 
numbers. The integral-method approach of references 7 and 8 was 
therefore studied and the results of the investigation conducted at 
the NACA Lewis laboratory are presented herein. 

The method of references 7 and 8 is based on an integral form of 
the equation of continuity .and on the intrinsic form in terms of the 
local streamline curvature of the irrotationality condition. The 
essential approximation of the method is the choice of a curvature 
function that will adequately represent the variation of the stream- 
line curvature in the flow field. In the two-dimensional problem, 
several guides are provided to limit the possible choices of a 
curvature function. These guides are; (a) the incompressible flow 
pattern; (b) the Prandtl-Glauert rule for the low -subsonic-speed 
range; and (c) the transonic similarity rule for the high-subsonic- 
speed range in which local supersonic regions exist in the flow field 
(references 9 and 10) . For the axially symmetric case, guides (a) and 
(b) are available but the situation with regard to guide (c) is dubious. 

An analysis of the differential equation for the body of revolution 
(reference 11) indicates that, although the transonic similarity rule 
of reference 9 may hold in the flow field, it does not appear to be 
valid on the 'body, where a different similarity rule obtains that is, 
moreover, valid only for extremely slender bodies . In this report, 
therefore, guide (c) has not been used in choosing a. curvature function. 
The resulting greater arbitrariness of the choice of curvature function 
thus implies a greater uncertainty of results for bodies of revolution 
than for the two-dimensional shapes. In order to reduce this uncertainty 
two substantially different curvature functions, both satisfying guides 
(a) and (b), have been selected for the investigation. 

Although the illustrative examples given herein have been confined 
to ellipsoids of revolution, the method with either curvature function 
may .be .applied in a straightforward manner to arbitrary bodies of 
revolution. 
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ANALYSIS 

The flow is assumed steady, isentropic, irrotational, and axially 
symmetric. Except for a modification in the continuity integral, the 
equations of motion and their development parallel the two-dimensional 
treatment of reference 8. The equations of motion are: (Main symbols 

are defined in the appendix; see fig. 1 for flow notation.) 

Equation of continuity, 


TO 


Pn 


Sitypo^o dy = 

JO. v 


2rtypv dn 


Irrotationality condition, 


dv 

5£ + 


Cv 


0 


( 1 ) 


( 2 ) 


Bernoulli equation, 


_ 1 



where 

p density of fluid 

v velocity in flow field 

n distance measured along potential lines from body of revolution 
C streamline curvature 
7 ratio of specific heats 
M Mach number 


subscript 0 denotes free-stream conditions. 

X % 

It is. extremely simplifying to replace the coordinate n by the 
coordinate y in the preceding equations. Because of the dubious 
situation with regard to a transonic similarity rule, it does not appear 
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easy to give a rigorous justification- for this replacement in the small- 
perturbation case analogous to the justification for the two-dimensional 
case given in appendix A of reference 8. The replacement of n by y, 
however, is exactly valid at the midchord location of bodies of revo- 
lution with fore-and-aft symmetry. The approximation involved in the 
replacement should therefore be fairly good near the maximum- thickness 
locations of slender bodies of revolution in general. For convenience, 
the free -stream quantities p Q and v Q will be taken as units of 

density and velocity so that p and v hereinafter represent non- 
dimensional fractions of free-stream density and velocity, respectively. 
In addition, the unit of length is taken as the airfoil chord, so that 
all lengths such as, x, y, and l/c hereinafter represent nondimens ional 
fractions of the airfoil chord length .- 

Equations (l), (2), and (3) thus become, respectively. 




°y 0 r y 

y dy =| pvy dy 

0 


(4) 


dy * "57 {5) 

1 

P = [l - Mo 2 (v 2 -l)] 7 1 (6) 

The variations indicated in equations (4) and (5) occur at constant x. 
The lower limit of integration Y in equation (4), which corresponds 
to the zero streamline, is the ordinate of the meridian section of the 
body of revolution. The variable upper limit y occuring in 
equation (4) and corresponding to the flow-field streamline may be 
eliminated by consideration of the limiting form of this equation as 


py 

y — > 00 . In order to obtain a finite 'result, the quantity y dy is 

Jy 

subtracted from both sides of equation (4). If no wakes or sources 
of fluid exist in the flow field, then lim (y-yo) = 0 (reference 8). 

, y — >oo 
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Equation (4) therefore assumes the limiting' form, for y — » 



r iOD 

' (pv-l) y dy 




Y 


(7) 


The continuity equation (7) becomes, when dy is replaced by dv 
from equation (5), 



( 8 ) 


in which the upper limit V is the velocity at the surface of the body 
corresponding to the ordinate Y. The relation between the flow -field 
velocity v and the corresponding distance y is obtained by inte- 
gration of equation (5) between the limits Y and y as 


y-Y = 


PV 


dv 

Cv 


O) 


For further simplification, the assumption of small perturbation - 
of the flow velocity is now explicitly introduced; accordingly, the 

quantity p - ^ in equation (8) is expressed as a power series in the 

* 1 
velocity increment o = v-1. From equation (6), the quantity p - — 

can be expanded in the form • 

P - ^ = UP - r c ^ 2 + r x u 3 (10) 

where 

U = 1 - Mq ( lla ) 

r c - 1 ■+ -§“ " ( m ) 

= 1 + ~21m 0 4 - | (2-y)(3-2r) Mq 6 (lie) 
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Terms of higher order than n 3 in the expansion equation (10) are 
negligible in the solution to be derived herein. Substituting 
equation (10) into the continuity integral equation (8) and changing 
to a new variable of integration z = u/ A yield 


— - A 2 

2 ~ A 


r i 


0 


(|iZ - r c Az 2 + r 1 A 2 z 3 ) £ dz 


( 12 ) 


where A= V-l is the disturbance velocity on the body. 

Equation (9) becomes, in terms of the new variable z, 

<^1 


y = Y + A 


Jz 


dz 

C(l+Az) 


(13) 


Inasmuch as 0 < z < 1 and for small perturbations A< < 1, the 
quantity _ may be expanded in powers of Az, yielding 


1+Az 


y = Y + A 


1 

U 


(1 - Az +A 2 z 2 - A 3 z 3 + • . .)dz 


(14) 


By inserting the value of y from equation (14) in equation (12), the 
continuity equation becomes', for arbitrary curvature C, 


Y2 

2 



^1 

/ ’ ' o ,\r P 1 i 

A 2 

i 

hiz-r c Az 2 +1i 2 JJY+/U i (l-Az+A 2 z 2 -A 3 z 3 + . . . )dz > 


0 ' 

\ L X 0 Jl 


I ' 


dz 

C 

'(15) 


First curvature function. - The streamline curvature is now 
assumed to be expressible as some simple function of the chosen inde- 
pendent variables, the chordwise location x, and the local velocity v. 
The function chosen must satisfy the boundary conditions of zero curva- 
ture at infinity and known curvature C a on the body. Two curvature 
functions were chosen for investigation. The first curvature function 
is the same as that used for the two-dimensional problem (reference 8). 



( 16 ) 
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where 

C a curvature of body contour 

p parameter to he determined 

The parameter p has a limited range of values because of conditions 
at infinity. The streamline curvature C must vanish at infinity and 
the integrals of equations (14) and (15) must diverge and converge, 

respectively, at the lower limit z = 0. Examination of equations (14) 

to (16) therefore shows that the parameter - p must have the range 
of values (See also reference 7.) 

1 < p c 3/2 (17) 


Substitution of equation (16) into equation (14) yields 


y = y + ~ 

u a 


l-z 1 ^ _ A (l-z 2 ' p ) A 2 (1 -z 5 -P) _ a 3 (1 -z 4 ~P) 
1-p 2-p 3-p ~ 4-p 


( 18 ) 


The continuity integral equation (12) becomes, by use of equation (18), 

n 

r z _ _ tij 
2 ~ c„ 


(^z 1 -? - IV z 2 "P + J y + A -A + 

v 1 c a L 1_ P 2 -P 


Uo 


A 


2 (l-z 3 ~ p ) a 3 (1 -z 4 ~ p ) 
3-p " 4-p 




>dz 


(19) 


Performing the integration and retaining terms through A 4 result in 


JjL. - Y0. + _J4 


2 A' 


,2-p 3-p 4-p 


(2-p)(3-2p) V (2-p)(4-2p) 


■jo 

(3-p)(4-2p| / 


( 20 ) 
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Several terms in equation (20) are negligible throughout the 
entire subsonic range. Thus, of the terms involving YC a as a factor, 

the term in is negligible when compared with the preceding one for 

A << 1. Of the terms not involving YC a as a factor, the quantity in 

p 

A is negligible when compared with the term in p A in the subsonic 
range p ~ 1, and in the A^ terms, the term in p is negligible 
compared with the term in F q in the transonic range p ~ 0. Hence, 

the p term within the second parenthesis on the right may be neglected 
and the continuity equation (20) becomes 




YC a + 


mA 

(2-p)(3-2p) 


t 0 a 2 

(3-p)(4-2p) 


( 21 ) 


Equation (21) provides the desired relation from which the velocity 
increment A may be obtained for given values of the free -stream Mach 
number Mq and the curvature parameter of the body YC & . Curves of A 

against ^ YC a for various values of Mq and the parameter p are 

presented in figure 2 . The numerical data from which these curves 
were drawn are presented in table I. 

i 

In applying equation (21), it is necessary first to obtain the 
parameter p. The parameter p is determined by requiring equation (21) 
to yield, in accordance with the discussion in the Introduction; (a) the 
known velocity A^ for incompressible flow p = 1, and (b) the Prandtl- 
Glauert rule for -bodies of revolution in the low-subsonic-speed range 
(p ~ 1). For incompressible flow, the terms in F Q may be neglected 

and equation (21) becomes, with p = 1 and A = A i , 


(YCa ) 2 _ YCa + A i 

g : 2 “ 2-p (2-p) (3-2p) 


( 22 ) ' 


The equation for determining p from the incompressible velocity A^ 
and the body-curvature parameter " YC a is then, from equation (22) 




(l+Ai) 

YC a 



(23) 


/ 
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The minus sign is chosen for the square root because the value of p 
may not exceed 3/2. In the low subsonic or Prandtl-Glauert range, the 
terms in r, are still negligible and equation (21) has the fom 


(TC a ) 2 _ nYC a ^ mA 
2 A 2 2 -P (2-p)(3-2p) 


(24) 


The quantity YC a is of order t 2 } where T is the thickness ratio. 

By multiplying the numerator and the denominator of the left side by p 2 , 
equation (24) may be expressed as 


\xA « f(pT 2 ,p) 


(25) 


The Prandtl-Glauert rule for bodies of revolution ' is (reference 1), 


pA = g(pT 2 ,x) 


(26) 


where the f motion g is defined by the incompressible solution 

= g(T 2 ,x) (27) 

In order that equation (25) agree with the Prandtl-Glauert rule 
(equation (26)), the parameter p must be chosen as a function of pT 2 
and x such that, for the incompressible case p = 1, this function 
reduces to that given in equation (23). Thus, in equation (23), 7 2 

is to be replaced by pT 2 in determining p. The parameter p is 
thus a function of Mach number for bodies of revolution, in contrast to 
the two-dimensional case where it is constant. When the value of p 
is obtained from low-speed considerations, the velocity increment A 
at high subsonic speeds is determined from equation (21). The detailed 
method of application of equation (21) to the calculation of velocities 
at high subsonic Mach numbers will be discussed in the section 
APPLICATION TO ELLIPSOID OF REVOLUTION . 

When the velocity increment A is determined, the local velocity 
u at any distance y in the flow field may be obtained from 
equation (18), which becomes,, on retention of terms through A 2 , 
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y - y 


1 A_ 
+ P-1 C a 



l A? 

2-P C a 


(l-z 2 -P) 


• (28) 


The variation of streamline curvature C with distance y may then he 
obtained from equations (16) and (28). 

Second curvature function. - As previously discussed, it was con- 
sidered desirable to repeat the analysis with a second curvature function 
substantially different from the first. In this manner, the reliability 
of the solution could be assessed more accurately. For greater simplicity 
the second curvature function was also chosen to. eliminate the factor y 
in the integrand of the continuity integral equation (12). The second 
curvature function was therefore chosen as 




(29) 


where the parameter q is to be determined as a function of x from 
the same low -speed considerations used to determine the parameter p 
in the preceding section. By an analysis similar to that used to 
determine the permissible values of p, the range of values of q is 
(reference 7), 

1 q < 2 (30) 

Substituting equation (29) into the continuity equation (12), integrating, 
and retaining terms through . yield 


2 • I= 2-q 3-q 


(31) 


Equation (31) provides the basic relation from which the velocity 
increment A may be calculated for given values of Mq and YC a . For 

incompressible flow, the term in F c is negligible so that the param- 
eter q is determined from 



2 2-q 


( 32 ) 
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In the Prandtl -Glauert range, the term in T c i 3 a tiH negligible, 
so that equation (31) becomes ^ ’ 


Yc a _ \iP? 
2 " 2-q 


(33) 


In order that equation (33) agree with the Prandtl -Glauert rule for 
bodies of revolution (equation (26)) the parameter q is chosen as a 
function of and x such that for the incompressible case this 

function reduces to equation (32). The velocity increment A is then 
determined at high subsonic speeds from equation (31), which may be 
conveniently solved by introducing two new parameters A and x (ref- 
erence 9). When A and X are defined by 



r / (2-q) 3 
2 m 3 (3- 9 ) 2 


(34a) 

(34b) 


equation (31) assumes the form 


X = A 2 - X 3 (35) 

A plot of equation (35) is presented in figure 3; the significance of 
the vertical tangents in figures 2 and 3 and their relation to the 
so-called potential -limit phenomenon are discussed in references 7 anrj 8. 

After the velocity increment A on the body is obtained, the local 
velocity increment v at any distance y in the flow field can be 
obtained by means of equation (5). Introducing equation (29) into 
equation (5) and changing over to the variable z yield 


py 


pi 


°a 


dz 

z*l(l4Az) 


(361 


Expanding the quantity l/l+Az, integrating term by term, and retaining 
terms through A^ give 
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y 2 _ y2 = 


2YA 1 ( 1 

C a q.-i \zl-l 


0 




(37) 


The variation of streamline curvature with distance y may then be 
obtained from equations (29) and (37). 


APPLICATION TO ELLIPSOID OF REVOLUTION 


As an illustration of the foregoing analysis, the velocity dis- 
tribution on an ellisoid of revolution at high subsonic Mach numbers 
will be calculated. The incompressible velocity distribution for an 
ellisoid of revolution of thickness ratio T may be obtained from the 
expression for the velocity potential given in reference 12 and is 
given by 


1 + A ± 


[l + f(T) 



l-£ 2 

| 2 (1-t 2 ) 


(38) 


where 


f(T) 


T 2 log© 

1 - V 1 -T 

2 a/ l-T 2 - T 2 logg ?- . + V A - lZ -T 

1- a/i-t 2 


(39) 


£ chordwise location measured from center of ellipsoid, semichords. 
The curvature parameter YCa for an ellipse is given by 


YC a = 


[l-e 2 (l-T 2 )] 3/2 


(40) 


First curvature function. - The parameter p is determined from 
equation (23) using the values of A± and YC a from equations (38) 
and (40), respectively. Although not so determined herein, the value 
of p may also be obtained directly from the curves of figure 2 for 
Mq = 0. In accordance with the discussion on the Prandtl-Glauert rule, 

p is to be determined as a function of pT 2 , which reduces to the 
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known dependence of p on t 2 for incompressible flow (p = l). 

The variation of p with p is conveniently obtained by plotting 
p against t 2 a nd changing the abscissa designation to pT 2 . The 
variation of p with p i 2 for several chordwise locations is pre- 
sented in figure 4. The values of p increase with increasing dis- 
tance from the center of the body and approach a limiting value of 
3/2. The occurrence of this limiting value may be seen by noting that 
the incompressible velocity passes through zero near the trailing 

edge and, from equation (23), the parameter ~p approaches 3/2 as 
approaches zero. Having obtained the value of p at the Mach number 
Mq, the velocity increment A for a given YC a is then obtained from 
the basic graphs of figure 2. Except for the additional step of finding 
p as a function of Mach number, the procedure is essentially the same 
as that given in reference 7. 


The parameter p has a limiting value of 3/2 for either Af — >0 
or p — >0. The graphical procedure for obtaining p therefore tends 
to become somewhat inaccurate at high Mach numbers unless great care 
is used, especially at those locations off the center where the velocity 
is close to zero. The value of p was therefore analytically deter- 
mined at the higher Mach numbers . Because the value of p at a Mach 
number Mo is the same as the incompressible value of p for a thick- 
ness ratio reduced by the factor */p, only higher -order terms were 
retained. Inasmuch as for p — »3/2, 2-p— *£■, equation (22) becomes, 

on retention of higher-order terms. 


3-2p 


4Aj 5 

(YC a ) 2 


(41) 


Equation (38) becomes, on retaining terms in 


f (t ) through 


t 2 


> 


1 + Aj. = 



l-l 2 


i4 2 (i-t 2 ) 


(42) 


Substituting equations (40) and (42) in equation (41) and replacing 

O O 

T by pT yields, for p as a f miction of Mach number, 

3 [l-fc 2 (l-MT 2 )] 3 

P 2 t 3 4 (1-£ 2 ) 

(43) 
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The velocity increment A was then f ound in the usual manner with the 
aid of the basic graphs (fig. 2). The procedure involving the analyti- 
cal determination of p was used in the range Mq >0.95. 


Second curvature function . - The parameter q is determined by 
means of equation (32) using the values of and YCa from 

equations (38) and (40), respectively. The variation of q with Mach 
number is then obtained by plotting q against t 2 and changing the 
abscissa designation to pT 2 . The variation of q with pT 2 for 
several chordwise locations is also presented in figure 4; the values 
of q are closer to the limiting value of 2 for the larger chordwise 
locations. When the value of q at the Mach number Mq is obtained, 
the parameter X is then calculated from equation (34b) . The value 
of the parameter k is next obtained from the X, X curve of figure 3. 
The velocity increment A .is finally obtained from equation (34a). 
Except for the additional step of finding q as a function of Mach 
number, the procedure is essentially the same as that of reference 8. 

It may be noted that calculations with the second curvature function 
are simpler than those with the first curvature function . 


As for the first curvature function, it is preferable to determine 
analytically the parameter q at high Mach numbers . Substituting 
equations (40) and (42) in equation (32) and replacing T 2 by pT 2 
yields, for q as a function of Mach number, 



,pT‘ 


-2\ I— 2dL 

)\ l \|i-$ 2 (i 


( 44 ) 

(1-pT 2 ) f pT 2 4l-£ 


The velocity increment A is then obtained by. the procedure following 
the graphical determination of q. The procedure involving the ana- 
lytical determination of q was used for Mach numbers Mq > 0.95. 

RESULTS AND DISCUSSION 

The methods developed herein for both curvature functions have 
been used to calculate velocity distributions on ellipsoids of revo- 
lution and conditions in the flow field for several thickness ratios 
and a range of Mach numbers. The variation of the maximum velocity 
increment Amax at the center of the ellipsoid with free -stream Mach 
number Mo for several thickness ratios is presented in figure 5. In 
order to indicate the local Mach number M corresponding to a given 
velocity increment A , contours of local Mach number obtained from 
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the Bernoulli relation between velocity increment A, local Mach 
number M, and free -stream Mach number Mo have been plotted on 
the figure. For the purpose of comparison, the values given in ref- 
erence 6 and the values given by the Prandtl-Glauert rule (equations - 
(26), (27), and (38)) are also included. For purposes of comparison, 
the Prandtl-Glauert rule is shown in figure 5 for local Mach numbers 
greater than 1, although this procedure is, strictly speaking, not 
permissible. On the whole, good agreement among the various methods 
is indicated. 

The velocity distribution on the body for several thickness ratios 
and Mach numbers is presented in figure 6. These distributions have 
been terminated at the values of t for 'which A passes through zero 
because the curvature functions chosen for this investigation, inasmuch 
as they yield positive values of the streamline curvature, may not be 
used for negative velocity increments. The agreement between the two 
curvature functions in figure 6 is essentially the same as that indi- 
cated in figure 5. The Mq i curves in figure 6 will be discussed 
in the section POTENTIAL LH4IT. 

The results with regard to the flow field are presented in fig- 
ures 7 to 9. The ' variation of the velocity parameter z with lateral 
distance y-T is given in figure 7. The curvature C/C a is plotted 
as a function of y-T in figure 8 and as a function of z in fig- 
ure 9. In contrast to the fairly good agreement for the velocity dis- 
tributions on the body, the results in the flow field are not in good 
agreement. This behavior is plausible inasmuch as the procedure for 
determining the curvature function utilizes the velocities on the body 
in th^ Prandtl-Glauert speed range as reference values. Thus, whereas 
the' method used herein to determine the curvature function can be 
regarded as averaging out local differences in the areas under the 
C/c a against y-T curves in figure 8 to yield about the same surface 
velocities (in accordance with equation (5)), the local differences in 
the flow field may be substantial. 

It is shown in figures 8 and 9 that the two curvature functions 
are quite different in type, as is evidenced by the values of C/C a 
greater than 1 exhibited by the second curvature function at high Mach 
numbers and low thickness ratios. Inasmuch as the curvature C/C a 
should probably not exceed 1, this phenomenon may be regarded as a 
criterion for eliminating this type of curvature function from consider- 
ation. The results obtained by the first curvature function are there- 
fore probably more reliable than those obtained by the second curvature 
function. . 
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The lateral and longitudinal extents of the symmetric -type local 
supersonic regions that develop at the higher Mach numbers, y s -T and 
| 3 , are presented as functions of free-3tream Mach number in figures 10 
and. 11, respectively. These values were obtained from the velocities 
in the flow field (fig. 7), the velocity distribution on the body 
(fig. 6), and the local M = 1 contour of figure 5. The two curvature 
functions are not in good agreement with regard to lateral extent of 
the local supersonic region, the discrepancy becoming more evident with 
decreasing thickness ratio. The two curvature functions are in much 
better agreement with respect to chordwise extent of the supersonic 
region, a result to be anticipated in view of the good agreement of the 
velocity distributions on the body (fig. 6). In order to indicate the 
differences between two-dimensional and axially symmetric flow with 
regard to extent of the supersonic region, the results for the elliptic 
cylinder, calculated by the method of reference 8 are included in 
figures 10 and 11. Although the supersonic region starts at a much 
higher Mach number for the ellipsoid of revolution than for. the elliptic 
cylinder, the lateral extent y s -T attains much higher values for the 
ellipsoid of revolution than for the elliptic cylinder. 


POTENTIAL LIMIT 

The calculation of two-dimensional compressible flows by the 
integral method (references 7 and 8) yields a free-stream Mach number, 
called the potential -limit Mach number Mo, I, beyond which, for a 
given thickness ratio, a continuous solution for the whole body cannot 
be obtained. The same phenomenon is found to occur in the treatment of 
bodies of revolution by the present method. The potential-limit Mach 
number corresponds to the points of vertical tangency of the A against 
i\j YCa curves of figure 2 or of the \ against X curve of figure 3. 
The value of M 0 _ ^ for the curvature functions investigated herein may 
therefore be obtained by application of the conditions 



«= 0 
d\ 


to equations (21) and (35), respectively. For the first curvature 


( 46 ) 
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function there is obtained. 


YC a 


3p 

A V 5-2p 
2-p f 2\i 
V- P 




(47) 


and for the second curvature function the result is 



For the ellipsoid of revolution and, probably more generally, for 
aerodynamic shapes, the potential limit occurs at the midchord section. 
(See reference 7 for discussion of two-dimensional case.) The potential- 
limit values of |a and A for the first curvature function were 
obtained by simultaneous solution of equations (21) and (47) at the 
mid chord location. Because of the variation of the parameter p with 
Mach number and the complicated form. of the continuity equation (21), 
the work was considerably more difficult than for the corresponding two- 
dimensional problem. It was found convenient to adopt an iteration 
procedure in which the values of n were obtained from the continuity 
equation (21) and the corresponding values of A from the potential- 
limit equation (47). The parameter . p was determined from the simpli- 
fied equation (43) evaluated at the midchord location, \ = 0. Initial 
values of |_i and A for the iteration were obtained as follows: 

At the potential limit, the value of |i is small so that the 
parameter p is close to its limiting values of 3/2. The continuity 
equation (21) and the potential-limit equation (47) may therefore be 
approximated by, respectively, 

( 49 ) 


3% ' T ( YOa + t A ) 


(50) 
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Equating the right sides of equations (49) and (50) yields 

1 1 

2 4 

A ■ (I fe) ^ (51) 

The value of T c at the potential limit is close to its limiting 
value of 1.2 for Mq = 1. Taking P c = 1.2, the value of A j may 
then he obtained from equation (51) and the value of pj, from either 
equation (49) or (50). The initial values of p 2 ^ l were 

within a few percent of the final iterated values. After the potential - 
limit Mach number Mo , I is obtained, the velocity distribution on the 
body and in the flow field can be calculated in the usual manner. 


. For the second curvature function, the potential -limit values of 
p and A may be more directly obtained. Insertion of equation (48) 
into equation (35) yields 


YC a To 2 (2-q) 3 4 

*- ( 3 -q)2 = 27 


(52) 


where the parameter q is given by the simplified equation (44) eval- 
uated at £ = 0. The quantity p, was then obtained from a plot of 
the parameter X against p and 1 the corresponding value of Aj, from 
equation (48). 

The velocity distribution on the body at the potential limit is 
given in figure 6. The velocity gradient at the midchord exhibits, as 
in the case of two-dimensional flow, a finite discontinuity (reference 7 
The results for the flow field at the potential limit are given in 
figures 12 to 14. The discrepancies previously noted between the two 
curvature functions become even more evident at the potential limit. 

The extreme maximum values of C/C a exhibited by the second curvature 
function indicate that it is probably very inaccurate near the potential 
limit. The potential-limit values of the lateral and longitudinal 
extents of the local supersonic region are indicated by circles in 
figures 10 and 11. (Some potential limit values of ya -t are outside 
range of fig. 10.) 

The variation of the potential -limit values of free -stream Mach 
number and maximum velocity increment with thickness ratio is given 
in figure 15 and the lateral and. longitudinal extents of the super- 
sonic region at the potential limit as a function of thickness ratio 
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-in figure 16. For comparison, the corresponding curves for the ellip- 
tic cylinder are also presented in' figures 15 and 16. The potential- 
limit Mach number Mq j is much higher for the ellipsoid of revolu- 
tion than for the elliptic cylinder and is very close to 1. These high 
values of Mq^j may be interpreted as follows: A body of revolution 

produces a smaller relative decrease of cross-sectional area between 
the body streamline A-A and an arbitrary streamline B.-B (fig. 1) 
than the corresponding two-dimensional shape. The consequent smaller 
velocity increments for the body of revolution result in higher values 
of the free -stream Mach number at which the supersonic region has 
started to develop (lower critical Mhch number).' Thus, it is plausible 
that the potential limit will occur at a much higher free -stream Mach 
number for the body of revolution than for the corresponding two- 
dimensional shape. It is surprising, however, that the values of Mo I 
are as close to 1 for a body of revolution as is indicated by the 
present method. The values of Mo j for the second curvature function 
are higher than those for the first curvature function. The difference 
in Mo, i for the two curvature functions may be explained on the basis 
of the large values of curvature exhibited by the second curvature 
function (figs. 9 and 13). These large values of curvature lead to a 
rapid decrease of velocity with distance from the body (equation (5)). 

It is therefore to be expected that the lateral extent of the local' 
supersonic region at a given Mach number would be smaller for the 
second curvature function than for the first curvature function. 

(Fig. 10 verifies this conclusion.) It then follows that the flow 
described by the second curvature function will not attain the maxi mum . 
value of Mo beyond which a symmetric -type solution may not be con- 
tinued (reference 7) until the flow given by the first curvature 
function has done so. 

The variation of maximum velocity increment at the potential limit 
i Wv , 7. with thickness ratio T for the first curvature function is 
approximately linear and close to that given by the initial choice of 
A may , i for iteration previously discussed (equation (51)). The vari- 
ation of A mo Y , i with t for the second curvature function is of the 
2/3 power type, as may.be seen from equation (48) and (52). This vari- 
ation, however, is the same as that given by the similarity law for two- 
dimensional flow (references 8 and 9), and illustrated by the curve of 
A max 2 against T for the elliptic cylinder. The description of 
A max , 7. hy the second curvature function is therefore to be regarded 
as another of its inaccuracies . 

The lateral extent of the supersonic region at the potential limit 
is much greater for the ellipsoid of revolution than for the elliptic 
cylinder. The extreme extent of the supersonic region for the ellipsoid 
of revolution is largely due to the closeness of Mo 2 to 1. For 
values of free-stream Mach number near 1, the disturbance velocity 
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decays very slowly with lateral distance from the "body (reference 13) . 
In addition, a very small velocity increment is needed to yield sonic 
velocity. The lateral extent of the supersonic region should there- 
fore rise very rapidly in the neighborhood of the potential limit . 

(See figs. 10 and 16.) The extreme behavior of the flow field in the 
neighborhood of the potential limit indicated by the present method 
may therefore be reasonable for potential flow. It is, of course, 
possible that the symmetric -type flow given by the present method 
would be terminated by a shock before such extreme sizes of the super- 
sonic region are attained. 

f 

On the basis of steady -state considerations, the potential-limit 
Mach number Mo j is interpreted in reference 8 as the transition Mach 
number between symmetric -type solutions for Mq ^ Mq^j and asymmetric- 
type solutions for Mo > Mo , j, with a shock terminating the symmetric- 
type solution. The extreme closeness of Mo , i to 1 would therefore 
indicate that the drag -divergence Mach number associated with the 
strong shock that terminates the symmetric -type local supersonic, 
region should occur much closer to 1 for bodies of revolution than for 
the corresponding. two-dimensional shapes. 


CONCLUDING REMARKS 

An integral method has been extended to obtain the compressible 
potential flow past a body of revolution. Of the two curvature func- 
tions chosen for the investigation, the first curvature function appears 
more reliable than the second curvature function. 

Although the symmetric -type local supersonic region first appears 
at a much higher Mach number for the ellipsoid of revolution than for 
the elliptic cylinder, it attains a much greater lateral extent. The 
potential -limit Mach number is much higher for the ellipsoid of revolu- 
tion than for the elliptic cylinder and very close to 1. The Mach 
number for drag divergence should therefore occur at a Mach number much 
closer to 1 for a body of revolution than for the corresponding two- 
dimensional shape-. The present approach offers fairly simple methods 
for calculating the velocity distribution on a body of revolution at 
high subsonic Mach numbers with local supersonic regions . 

Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 

Cleveland, Ohio, September 1, 1950. 
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APPENDIX - SYMBOLS 

The following symbols are used in this report : 

C streamline curvature 

M Mach number 

n distance measured along potential lines from body of revolution 
p,<l par ame ters for curvature functions 

V resultant velocity on body of revolution 
v velocity in flow field , 

x distance along axis of body of revolution measured from leading 
edge 

Y ordinate of meridian section of body of revolution 

y lateral distance normal to free-stream direction 

z = u/ A 



y ratio of specific heats 

A velocity increment, V-l 



p = 1 -Mq 2 


V velocity increment, v-l 

i distance along axis of body of revolution measured from center, 
s emi chord 

p density of fluid 

t thickness ratio 

2 I-* 3 (3-<l) 2 


X 
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Subscripts: 

a at surface of body of revolution 

i incompressible or low-speed value 

l potential limit 

s sonic value 

0 free-stream condition 
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TABLE I - BASIC CALCULATIONS FOR OBTAINING VELOCITY INCREMENT A (EQUATION (21)) 
[First curvature function, C/C a * 2 ?.]] 




A 

0.02 

0.04 

0.06 

0.08 

0. 10 

0.12 

0.16 

0.20 

0.25 

0.30 

0.40 

0.50 

0.60 

0.70 

0.80 

0.90 

i^ol 

*0 

P 

V YC » : 1 

0 

1.0 

1.1 

1.2 

1.3 

1.4 
1.45 
1.48 
1.495 
1.499 
1.4999 

0.0665 

.0718 

.0792 

.0895 

.1099 

.1344 

.1688 

.2389 

.3567 

.6331 

0.1138 

.1230 

.1355 

.1538 

.1881 

.2260 

.2855 

.4029 

.6005 

1.0652 

0.1563 

.1689 

.1859 

.2109 

.2570 

.3086 

.3886 

.5473 

.8146 

0.1960 

.2118 

.2329 

.2639 

.3210 

.3848 

.4838 

.6802 

1.0116 

0.2338 

.2526 

.2776 

.3144 

.3818 

.4569 

.5736 

.8053 

0.2702 

.2918 

.3206 

.3628 

.4399 

.5259 

.6593 

.9247 

0.3396 

.3666 

.4025 

.4550 

.5506 

.6568 

:8217 

1.1499 

0.4055 

.4377 

.4804 

.5426 

.6555 

.7806 

.9749 

0.4841 

.5225 

.5732 

.6471 

.7807 

.9280 

1.1569 

0.5592 

.6037 

.6623 

.7474 

.9004 

1.0691 

0.7006 

.7567 

.8305 

.9370 

1.1276 

1.3363 

0.8315 

.8991 

.9877 

1.1148 

0.9530 

1.0319 

1.1352 

1.2830 

1.0656 

1.1558 

1.2735 

1.1693 



0.4 

1.0 

1.1 

1.2 

1.3 

1.4 
1.45 
1.48 
1.495 
1.499 
1.4999 

0.0628 

.0688 

.0756 

.0852 

.1048 

.1284 

.1614 

.2286 

.3413 

.6061 

0.1082 

.1169 

.1290 

.1466 

.1795 

.2158 

.2729 

.3854 

.5745 

1.0194 

0.1483 

-.1604 

.1766 

.2007 

.2449 

.2944 

.3712 

.5234 

.7795 

0.1857 

.2009 

.2210 

.2508 

.3057 

.3669 

.4621 

.6505 

.9680 

0.2213 

.2392 

.2632 

.2985 

.3633 

.4354 

.5476 

.7699 

0.2553 

.2760 

.3036 

.3441 

.4183 

.5009 

.6292 

.8838 

0.3200 

.3459 

.3803 

.4308 

.5228 

.6250 

.7837 

1.0988 

0.3811 
.4120 
.4529 
.5128 
.6217 
. 7423 
.9392 

0.4534 

.4903 

.5393 

.6102 

.7392 

.8816 

0.5221 

.5650 

.6213 

.7033 

.8514 

1.0144 

0.6498 

.7038 

.7752 

.8782 

1.0630 

0.7660 

.8313 

.9170 

1.0406 

0.8715 

.9479 

1.0482 

0.9666 

1.0545 

1.1694 

1.0515 



0.6 

1.0 

1.1 

1.2 

1.3 

1.4 
1.45 
1.48 
1.495 
1.499 
1.4999 

0.0584 

.0634 

.0700 

.0793 

.0976 

.1197 

.1506 

.2134 

.3187 

.5661 

0.1000 

.1083 

.1197 

.1359 

.1668 

.2009 

.2543 

.3596 

.5366 

.9523 

0.1369 

.1482 

.1634 

.1859 

.2273 

.2738 

.3458 

.4883 

.7278 

0.1709 

.1851 

.2040 

.2319 

.2834 

.3409 

.4301 

.6066 

.9037 

0.2032 

.2200 

.2424 

.2755 

.3364 

.4042 

.5182 

.7178 

0.2339 

.2533 

.2791 

.3171 

.3869 

.4647 

.5851 

.8238 

0.2919 

.3162 

.3485 

.3959 

.4827 

.5790 

.7282 

1.0238 

0.3462 

.3752 

.4137 

.4701 

.5729 

.6867 

.8628 

0.4099 

.4445 

.4904 

.5577 

.6796 

.8143 

1.0223 

0.4694 

.5095 

.5628 

.6405 

.7810 

.9356 

0.5776 

.6287 

.6963 

.7945 

.9707 

1.1632 

0.6725 

.7347 

.8168 

.9350 

1.1463 

1.3755 

0.7549 

.8284 

.9249 

1.0636 

1.3092 

1.5742 

0.8246 

.9096 

1.0217 

1.1811 

1.4615 

1.7618 

0.8798 

.9794 

1.1067 

1.2876 

1.6031 

1.9388 

0.9257 

1.0367 

1.1807 

1.3836 

1.7355 

2.0164 

0.9563 

1.0817 

1.2430 

1.4694 

1.8580 

2.2652 

0.7 

1.0 

1.1 

1.2 

1.3 

1.4 
1.45 
1.48 
1.495 
1.499 
1.4999 

0.0551 

.0596 

.0662 

.0747 

.0920 

.1129 

.1421 

.2015 

.3011 

.5349 

0.0939 

.1017 

.1122 

.1278 

.1570 

.1893 

.2399 

.3396 

.5068 

.8996 

0.1281 
.1387 
. 1531 
.1744 
.2137 
.2578 
.3260 
.4609 
.6873 
1.1948 

0.1596 

.1730 

.1909 

.2174 

.2662 

.3207 

.4053 

.5636 

.8532 

0.1893 

.2052 

.2265 

.2579 

.3156 

.3801 

.4799 

.6772 

1.0091 

0.2176 

.2359 

.2604 

.2965 

.3628 

.4367 

.5510 

.7771 

0.2705 

.2935 

.3241 

.3692 

.4518 

.5435 

.6853 

.9654 

0.3196 

.3471 

.3837 

.4374 

.5354 

.6439 

.8114 

0. 3764 
.4095 
.4533 
,5175 
.6340 
.7626 
.9605 

0.4290 

.4672 

.5183 

.5926 

.7271 

.8749 

0.5221 

.5713 

.6364 

.7308 

.9004 

1.0852 

0.6003 

.6607 

.7401 

.8548 

1.0589 

1.2796 

0.6642 

.7362 

.8305 

.9659 

1.2047 

1.4609 

0.7131 

.7978 

.9081 

1.0648 

1.3389 

1.6304 

0.7469 

.8455 

.9727 

1.1520 

1.4684 

1.7893 

0.7647 
.8789 
1.0245 
1.2277 
1.6754 
1 . 9383 

0.7658 

.8976 

1.0633 

1.2919 

1-.6785 

2.0779 

0.8 

1.0 

1.1 

1.2 

1.3 

1.4 
1.45 
1.48 
1.495 
1.499 
1.4999 

0.0509 

.0546 

.0603 

.0683 

.0842 

.1032 

.1300 

.1845 

.2759 

.4502 

0.0852 

.0920 

.1019 

.1162 

.1431 

.1728 

.2193 

.3108 

.4642 

.8244 

0.1157 

.1256 

.1387 

.1583 

.1945 

.2351 

.2978 

.4217 

.6296 

1.1178 

0.1437 

.1560 

.1725 

.1969 

.2419 

.2922 

.3708 

.5237 

.7815 

0.1699 

.1845 

.2041 

.2332 

.2865 

.3460 

.4381 

.6194 

.9241 

0.1946 
.2115 
.2341 
.2675 
.3289 ! 
.3972 
.5027 
.7106 

0.2404 

.2616 

.2901 

.3319 

.4086 

.4935 

.6245 

.8824 

0.2821 

.3076 

.3417 

.3917 

.4830 

.5837 

.7387 

1.0436 

0.3294 

.3601 

.4011 

.4611 

.5701 

.6899 

.8736 

0.3716 

.4078 

.4556 

.6256 

.6519 

.7901 

1.0011 

0.4424 

.4896 

.5516 

.6416 

.8021 

.9760 

0.4950 

.5541 

.6314 

.7421 

.9371 

0.5287 

.6016 

.6956 

.8281 

1.0587 

0.5419 

.6313 

.7440 

.9005 

0.5324 

.6424 

.7768 

.9915 

0.4965 

.6336 

.7935 

1.0047 

0.4287 

.6033 

.7940 

1.0371 

0.85 

1.0 

1.1 

1.2 

l;3 

1:4 

1:45 

1:48 

1.495 

1:499 

1.4999 

0.0459 

.0507 

.0560 

.0635 

.0786 

.0965 

.1217 

.1728 

.2584 

.4592 

0.0790 
;0855 
:0949 
:1083 
.1336 
: 1615 
:2051 
.2910 
:4347 
.7723 

0:1072 

:1164 

.1288 

.1474 

.1814 

.2195 

.2785 

.3947 

.5896 

1.0472 

0.1327 

.1444 

.1599 

.1829 

.2253 

.2727 

.3459 

.4801 

.7318 

0.1566 

.1704 

.1889 

.2162 

.2666 

.3227 

.4093 

.5797 

.8654 

0.1789 
. 1948 1 
.2101 
.2477 
.3057 
.3701 
.4695 
.6648 
.9924 

0:2196 

.2398 

.2667 

.3064 

.3790 

.4593 

.5828 

.8252 

0:2661 

.2804 

.3129 

.3605 

.4472 

.5426 

.6889 

.9757 

0.2963 

.3260 

.3652 

.4226 

.5266 

.6403 

.8139 

0:3311 
.3661 
.4124 
.4797 
.6006 
. 7321 ! 
.9319 

0.3844 

.4312 

.4923 

.5800 

.7349 

.9013 

0.4151 

.4759 

.5539 

.6635 

.8536 

1.0548 

0.4201 

.4995 

.5972 

.7313 

.9583 

1.1945 

0.3933 

.4993 

.6219 

.7859 

1.0500 

0.3206 

.4716 

.6269 

.8205 

0.1445 

.4101 

.6111 

.8423 

0.2892 

.2992 

.5731 

.8490 

0.9 

1.0 

1.1 

1.2 

1.3 

1.4 
1.45 
1.48 
1.495 
1:499 
1.4999 

0.0425 

.0456 

.0504 

.0572 

.0716 

.0876 

.1105 

.1571 

.2350 

.4179 

0.0712 

.0771 

.0857 

.0977 

.1208 

.1463 

.1862 

.2644 

.3955 

.7026 

0.0958 
. 1043 
.1158 
.1326 
.1638 
.1987 
.2526 
.3586 
.5361 
.9524 

0.1182 

.1288 

.1430 

.1642 

.2032 

.2466 

.3136 

.4451 

.6653 

0.1387 

.1513 

.1684 

.1937 

.2400 

.2915 

.3708 

:5263 

.7866 

0.1575 

.1723 

.1921 

.2213 

.2748 

.3341 

.4251 

.6034 

.9034 

0.1912 
.2101 
.2353 
.2722 
.3396 
.4138 
: 5272 
.7487 

0.2200 

.2430 

.2737 

.3184 

.3994 

.4878 

.6224 

.8848 

0.2497 
.2782 
. 3160 
.3705 
.4682 
.5785 
.7343 
1.0452 

0.2722 

.3069 

.3522 

.4170 

.5317 

.6591 

.8396 

0.2937 

.3443 

.4074 

.4947 

.6443 

.8072 

1.0347 

0.2732 

.3512 

.4388 

.5528 

.7404 

.9394 

0.1591 

.3170 

.4439 

.5916 

.0212 

1.0573 

0.2049 

.4176 

.6107 

.8875 

1.1626 

0.3494 

.6093 

.9401 

1.2561 

0.2047 

.5863 

.9791 

1.3113 

0.5310 

1.0053 

1.4097 
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TABLE I - BASIC CALCULATIONS FOR OBTAINING VELOCITY INCREMENT A (EQUATION - (21 ) ) - Concluded. 
[First curvature function, C/C fl = z p /] 



A 

0.02 

0.04 

0.06 

|- 0.08 

0.10 

0.12 

0.16 

! 0.20 

0.25 

0.30 

0.32 1 0.3218 

0.36 

0.30 

0.384 

*0 

P 



• w: — • 1 

0.95 

1.0 
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(a) Free-stream Mach number M Qf 0. 

Figure 2. - Velocity increment on surface of body of revolution as. function of curvature 
parameter of body V^a ** or several values of parameter p. Values computed from 
equation (21) and tabulated in table I. 
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Curvature parameter, VYC^ 


(b) Free-stream Mach number Mq, 0.4. 


Figure 2. - Continued'. Velocity increment on surface of body of revolution . as function of 
curvature parameter of body for several values of parameter . p. Values computed 

from equation (21) and tabulated in table I. 





NACA TN 2245 



Figure 2. - Continued. Velocity increment on surface of body of revolution as function of 
curvature parameter of body V YC a for several values of parameter p. Values computed 
from equation (21) and tabulated in table I. 
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Curvature parameter, VYC a 
(d) Free stream Mach number Mg, 0.7. 

Figure 2. - Continued. Velocity increment on surface of body of revolution as function of 
curvature parameter of body ^YC a for several values of parameter p. Values computed 
• from equation (21) and tabulated in table I. 
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Figure 2. - Continued. Velocity increment on surface of body of revolution as function of 
curvature parameter of body £C a for several values of parameter p. Values computed 
from equation (21) and tabulated in table I. 


Velocity increment on body 
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Figure 2. - Continued. Velocit y in crement on surface of body of revolution as function of 
curvature parameter of body C a for several values of parameter p. Values computed 
from equation (21) and tabulated in table I. 



Velocity Increment on body. 
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(g) Free-stream Mach number Mq, 0.90* 

Figure 2. - Continued. Velocit y In crement on surface of body of revolution as function of 
curvature parameter of body Y YC a for several values of parameter p. Values computed 
from equation (21) and tabulated In table I. 



•Velocity increment on body. 
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Curvature parameter, VYC a 

(h) Free-stream Mach number Mq, 0.95. 


Figure 2. - Continued. Velocit y inc rement on surface of body of revolution as function of 
curvature parameter of body aJxC ^ for several values of parameter p. Values computed 
from equation (21) and tabulated in table I. 




Velocity increment on body. 
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Figure 2. - Continued. Velocity increment on surface of body of revolution as function of 
curvature parameter of body for several values of parameter p. Values computed 

from equation (21) and tabulated in table I. 



Velocity increment on body. 
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Figure 2. - Concluded. Velocity increment^ on surface of body of revolution as function of 
curvature parameter of body* YYC a for several values of parameter p. Values computed 
from equation (21) and tabulated in table I. 
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Figure 3. - Plot of equation (35); X » 
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Figure 4. - Curvature-function parameters p and q as functions of p 
for several chordwise positions on ellipsoid of revolution. 



Maximum velocity increment 
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Figure 5. - Variation of maximum velocity Increment at midchord of ellipsoid of 
revolution as function of free-stream Mach number for several thickness ratios. 
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Figure 7. - Lateral variation of velocity in flow field of ellipsoid of revolution at midchord location for several 

Mach numbers. 
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Figure 8* - Lateral variation of streamline curvature in flow field of ellipsoid of revolution. 

at midchord location for several Mach numbers® 
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Figure 3. - Concluded. Lateral variation of streamline curvature in flow field of ellipsoid; 
of revolution at midchord location for several Mach numbers. 




Eatio of streamline curvature to body curvature, C/C 
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(b) Thickness ratio t, 0.3. 

Figure 9. - Concluded. Variation of streamline curvature with velocity in 
flow field of ellipsoid of revolution for several Mach numbers. 
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Figure 10. - Lateral extent of supersonic region at midchord location for ellipsoid 

of revolution and elliptic cylinder. • 






First curvature function 
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Figure 11. - Longitudinal extent of supersonic region for ellipsoid of revolution and elliptic cylinder, 







50 


NACA TN 2245 



2 ‘.za^at nared X^pootaA p-tapj-Aop^ 


Figure 12. - Lateral variation of velocity in flow field of ellipsoid of revolution at mid chord location for 

potential -limit Mach number and several thickness ratios. 




Figure 13. - Lateral variation of streamline curvature in flow field of ellipsoid of revolution at midchord location 

for potential- limit Mach number and several thickness ratios. 


Ratio of streamline curvature to body curvature, C/C 
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Figure 15. - Potential-limit values of free-stream Mach number and velocity increment 
at midchord for ellipsoid of revolution and elliptic cylinder as function of 
thickness ratio* 
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Longitudinal extent of supersonic region 
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Thickness ratio, t 

Figure 16. - Lateral and longitudinal extents of supersonic region at potential 
limit for ellipsoid of revolution and elliptic cylinder as function of thickness 
ratio. 
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